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I  INTRODUCTION 


The  objectives  of  this  project  are  to  develop  the  theoretical  tools 
necessary  to  treat  a  variety  of  molecular  state-to-state  collision  processes 
and  to  perform  calculations  of  transition  probabilities  for  selected 
systems.  The  problems  considered  here  have  involved  the  interactions  of  three 
particles  under  the  influence  of  potentials  for  two  or  more  electronic 
states.  For  this  type  of  process ,  currently  available  experimental  techniques 
allow  the  determination  of  specific  states  of  species  before  and  after  inter¬ 
action.  Our  general  approach  has  been  to  pursue  two  types  of  calculations: 
model  calculations  designed  to  Isolate  and  probe  particular  reaction 
mechanisms,  and  calculations  focused  on  specific  systems  that  have  been 
studied  experimentally  in  our  laboratory  or  elsewhere.  This  dual  approach  has 
enabled  us  to  better  understand  a  wide  range  of  scattering  phenomena  and  has 
provided  a  quantitative  test  of  the  accuracy  of  the  methods  we  have  developed 
and  of  frequently  used  approximations. 

During  the  past  year  we  have  devoted  the  major  part  of  our  effort  to 
detailed  time-independent  quantum  mechanical  calculations  of  ion  pair  forma¬ 
tion  and  related  vibrational  excitation  processes  in  atom-diatomic  colli¬ 
sions.  We  also  are  applying  the  methods  developed  in  this  work  to  the  case  of 
symmetric  molecular  charge  exchange.  In  related  work,  we  have  treated  the 
process  of  ion-pair  formation  in  nonadlabatic  collisions  by  quasi-classical 
trajectory  surface  hopping  and  quantum  mechanical  perturbation  theory  tech¬ 
niques.  Finally,  previous  rigorous  coupled  channel  calculations  of  cross 
sections  for  1-mixing  of  Na  Rydberg  states  in  collisions  with  atoms  and 
molecules  have  been  complemented  by  the  development  of  a  simple  scaling  law 
for  these  cross  sections.  All  of  these  accomplishments  are  described  briefly 
in  the  following  section. 


II  RESEARCH  ACCOMPLISHMENTS 

Coupled-Channel  Calcwilationa  of  Vibrational  Excitation  apd  Ion-Pair  Formation: 

K.  +  0, 

A  significant  fraction  of  our  effort  under  this  program  has  been  devoted 
to  the  investigation  of  electronic-vibrational-translational  energy  exchange 
in  nonadiabatic  collisions.  Fully  quantum  mechanical  coupled  channel 
scattering  calculations  were  performed  on  a  model  potential  surface  designed 
to  model  the  system 


K  +  °2(v1)  ♦  K  +  02(vf ) 


♦  K  +  02  (vf) 


This  system  is  an  attractive  candidate  for  theoretical  study  because  experi¬ 
mental  data  are  available  at  different  collision  energies  by  different 
12  3 

workers,  ’  and  data  from  this  laboratory  are  also  available  for  the  related 
Ar*  +  02  system.  The  latter  results  showed  the  interesting  effect  of  media¬ 
tion  of  the  vibrational  distribution  in  the  initial  channel  by  the  presence  of 
the  ion-pair  channel.  Since  the  present  quantum  calculations  are  essentially 
dynamically  exact,  they  will  serve  as  a  benchmark  for  other  congmtatlonal 
approaches.  They  also  help  elucidate  the  dynamics  of  the  efficient  flow  of 
energy  between  different  vibrational  levels  in  the  presence  of  nonadiabatic 
coupling.  Finally,  in  the  course  of  making  these  calculations,  we  have 
Improved  the  state-of-the-art  methodology  for  streamlining  large  scale  scat¬ 
tering  calculations. 

In  these  calculations,  only  the  doublet  potential  surfaces  arising  from 

the  K(^S)  +  0_(^I  )  ground  state  asynptote  and  the  K+(*S)  +  0-  (^n  )  ion 
-  8  »  8 
pair  asysptote  were  considered.  Simplified  but  realistic  interaction  poten- 

o 

tials  were  taken  from  a  classical  collision  model  fit  to  data.  The  interac¬ 
tion  was  taken  to  be  isotropic  (no  rotational  excitation)  to  make  the  problem 
more  tractable  and  to  focus  on  the  vlbronlc  aspects  of  the  problem.  Coulomb 
boundary  conditions  for  the  ion-pair  formation  channels  had  to  be  incorporated 
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explicitly.  The  collision  energies  studied  were  over  the  broad  range  0.27, 
0.68,  12.5,  and  43.5  eV.  Predominantly  a  basis  set  was  used  consisting  of 
16  neutral  and  16  ionic  vibrational  eigenfunctions.  Extensive  testing  of 
convergence  with  number  of  vibrational  channels  was  performed. 

The  results  of  the  calculations  show  that  vibrational  energy  transfer  is 
efficient  in  nonadiabatic  collisions,  typically  being  a  significant  fraction 
of  the  gas-kinetic  cross  section.  At  the  two  highest  energies,  the  integral 
cross  sections  (Q)  for  scattering  from  K  +  (^(v  -  0)  for  ion-pair  formation 
are  very  roughly  governed  by  the  Franck-Condon  factors  for  the  transi¬ 

tions.  The  neutral  pair  Q's  show  a  slow  monotonlc  fall-off  with  v.  At  low 
energies,  the  Q's  are  dominated  by  energy  availability.  The  effect  of  0£~ 
vibration  can  be  seen  in  the  opacity  functions  and  differentirl  cross  sections 
that  were  calculated  quantum  mechanically,  and  the  effect  can  be  related  to 
the  collision  time  and  the  0£~  vibrational  period  in  a  classical  surface¬ 
hopping  model.  The  structure  in  the  high  energy  loss  peaks  of  neutral 
products  from  TOF  measurements1  is  rather  well  reproduced  in  these  calcula¬ 
tions.  This  suggests  that  much  of  this  Inelasticity  is  due  to  ground  elec¬ 
tronic  state  vibrational  excitation. 

A  paper  on  this  work,  published  in  the  Journal  of  Chemical  Physics,  gives 
the  theoretical  methods  and  results  in  detail.  (See  list  of  papers  in  Section 
III  and  abstract  in  Appendix  A). 


>tric  Molecular  Charge  Exchar 


As  a  further  probe  of  efficient  vibrational  energy  exchange  in  nonadia¬ 
batic  collisions,  we  have  also  completed  a  study  of  symmetric  molecular 
scattering  in  the  ground  electronic  state  system 


02+(v1)  +  02(vp  -*■  °2+(yf)  +  °2^VP 

*  °2^VP  +  °2+^vP 


(2) 


The  direct  and  charge  transfer  (exchange)  scattering  characteristic  of  symmet¬ 
ric  systems  leads  to  much  larger  vlbratlonally  inelastic  and  momentum  transfer 
cross  sections  than  for  asymmetric  systems. 


One  previous  theoretical  formulation  for  this  problem  exists,  given  by 
Bates  and  Reld.^  It  approximates  the  treatment  of  translational  motion  by 
using  a  classical  trajectory.  Also,  to  make  the  problem  more  tractable,  the 
rotational  degree  of  freedom  was  ignored.  This  formulation^  is  a  form  of  the 
so-called  impact  parameter  approximation  of  scattering  theory.^  This  method 
is  frequently  applied  at  hyperthermal  and  even  thermal  energies,  often  without 
due  consideration  of  its  accuracy.  It  is  known  from  general  considerations 
that  the  approximation  improves  at  higher  energies,  and  it  is  usually  consid¬ 
ered  valid  for  most  cases  around  10  eV.^ 

A  completely  quantum  mechanical  formulation  for  symmetric  molecular  scat¬ 
tering  has,  to  our  knowledge,  been  lacking  until  now.  It  appeared  worthwhile 
to  examine  the  reliability  of  the  impact  parameter  approximation  (using  dif¬ 
ferent  classical  trajectories)  by  comparing  with  fully  quantum  mechanical 
coupled-channel  scattering  calculations.  Therefore,  we  formulated  the  problem 
and  performed  the  coupled-channel  calculations  at  center-of-mass  collision 
energies  of  1,  8,  and  36  eV. 

A  unitary  transformation  was  found  to  reduce  the  2n  coupled  radial  equa¬ 
tions  to  two  sets  of  n  coupled  equations.  This  reduction  in  dimensionality 
significantly  reduces  the  computational  labor  because  the  computing  time 
scales  approximately  as  the  cube  of  the  number  of  coupled  equations. 

The  coupled  radial  SchrtJdinger  equations  were  numerically  integrated  by 
the  log-derivative  algorithm  developed  by  Johnson.**’ ^  Convergence  was  checked 
for  Integration  step  size,  integration  range,  and  number  of  channels  needed  in 
the  basis  set.  It  was  found  that  basis  size  convergence  was  more  difficult  at 
1  eV  than  at  8  or  36  eV.  A  somewhat  similar  result  was  found  in  the  K.  +  O2 
coupled  channel  calculations. 

A 

The  impact  parameter  method  of  Bates  and  Raid  has  been  applied  by  Horan 
and  coworkers®-10  to  the  0^-02  system,  using  both  straight  line  and  effective 
curved  trajectories.  Comparison  of  the  relative  probabilities  for  transitions 
at  a  particular  impact  parameter  for  1  eV  collisions  from  the  present  calcula¬ 
tions  with  the  strai^it  line  and  curved  trajectory  classical  path  calculations 
of  Ref.  10  show  essentially  no  agreement  between  any  of  the  three  methods. 
Comparison  of  integral  cross  sections  shows  that  agreement  is  not  particularly 


close  between  any  of  the  approaches  at  8  eV  and  Is  qualitative  only  at 
36  eV.  It  Is  expected  that,  as  the  energy  Increases,  the  classical  path  ap¬ 
proximation  will  be  increasingly  accurate. 


A  paper  on  this  subject  is  to  be  published  in  the  Journal  of  Chemical 
Physics  (see  Appendix  B  and  list  of  publications  in  Section  III). 

Model  Calculations  of  Nonadiabatic  Collisions 

Under  this  program  we  have  also  completed  studies  of  atom-molecule 
collisions  that  involve  two  or  more  electronic  potential  energy  surfaces  using 
both  qua8i-classical  and  quantum  mechanical  perturbation  techniques,  as 
described  below. 


Quasi-Classical  Surface  Hopping  Model 

This  model  invokes  the  Landau-Zener  formula  to  calculate  a  surface¬ 
hopping  probability  at  each  surface  intersection.  Using  the  surface-hopping 
model,  we  performed  a  comparative  study  of  ion-pair  formation  on  two  similar 
systems: 


and 


Ar  +  I2 
K  +  I2 


-*•  Ar+  +  I2 

(3) 

+  K+  +  I2” 

(4) 

One  would  expect  these  two  reactions  to  be  rather  similar,  because  in  each 
case  the  single  active  electron  is  in  a  4s  orbital.  However,  an  additional 
final  channel  is  available  in  the  first  reaction.  Penning  ionization  may 
occur. 


Ar*  +  I2  Ar  +  I2+  +  e“  (5) 

11  12  13 

Data  are  available  for  both  ion-pair  reactions.  *  *  We  were  able  to 

perform  calculations  that  explained  some  of  the  differences  in  the  differen¬ 
tial  cross  sections  for  these  two  systems.  We  explored  how  differences  in  the 


matrix  elements  arising  from  the  inner-shell  vacancy  and  the  competing  ioniza¬ 
tion  channel  would  affect  the  cross  sections.  A  paper  on  this  work  has  been 
published  in  the  Journal  of  Chemical  Physics  (see  list  of  papers  in 
Section  111  and  abstract  in  Appendix  A). 

Quantum  Mechanical  Perturbation  Theory  of  Honadiabatic  Collisions 

In  this  work,  we  sought  an  alternative  to  the  coupled-channel  quantum 
mechanical  formulation  of  nonadiabatic  collisions  that  would  preserve  the 
essential  features  of  the  multiple-surface  crossings,  but  exploited 
approximtions  that  can  legitimately  be  made  at  high  collision  energies 
(E  ~  50-100  eV).  The  approach  chosen  was  the  impact  parameter  approximation, 
in  which  a  rectilinear  classical  trajectory  is  assumed  for  the  translational 
motion,  leading  to  a  time-dependent  SchrOdinger 's  equation  for  the  remaining 
degrees  of  freedom.  We  discovered  that  model  potentials  of  the  same  type  used 
in  previous  studies  of  collisions  of  alkali  atoms  and  diatomic  molecules  led 
to  a  particularly  simple  form  of  the  final  equations.  These  equations  were 
solved  using  the  Magnus  approximation. 

We  will  summarize  the  results  of  our  calculations.  We  performed  calcula¬ 
tions  using  model  potentials  chosen  to  model  the  reaction 

Cs  +  02  Cb+  +  0  "  (6) 

Initially,  isotropic  potentials  were  assumed;  that  is,  molecular  rotation  was 

neglected.  The  results  exhibited  strong  oscillations,  corresponding  qualita- 

14 

tively  to  those  observed  in  recent  measurements.  These  oscillations  have 
been  related  to  the  vibrational  motion  of  the  molecular  potentials  using  the 
infinite-order-sudden  approximation.  We  obtained  rotatlonally  averaged  cross 
sections  between  specific  vibrational  states.  The  oscillations  observed  in 
the  Isotropic  case  persisted,  although  they  were  less  prominent. 

These  rotatlonally  averaged  calculations  also  enabled  us  to  examine  the 
effect  of  different  initial  vibrational  levels  v^  on  the  ion-pair  formation. 
The  model  predicted  that  the  oscillations  will  become  less  pronounced  as  v^  is 
Increased.  Although  no  measurements  are  available  to  test  this  result,  our 
calculations  suggest  a  direction  for  future  experiments. 
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The  theory  and  calculations  are  described  in  a  paper  that  has  been  pub¬ 
lished  by  the  Journal  of  Chemical  Physics  (see  Section  III  and  Appendix  A) . 

Theoretical  Studies  of  Rydberg  Atoms 

During  the  past  year,  we  have  completed  a  study  of  physical  scaling  laws 
for  the  angular  momentum  mixing  of  Na  Rydberg  states  in  collisions  with  atoms 
and  molecules.  We  have  found  that  the  results  of  many  coupled-channel  and 
Born  calculations  can  be  expressed  in  terms  of  numerically  determined  func¬ 
tions  of  reduced  parameters.  The  general  result  is 

aA-tnix(n)  ‘  TO4ao2g(p)f(Y)  (7) 

where  a ^  ^  (n)  is  the  experimentally  measured  A—mixing  cross  section  summed 
over  all  final  states  A';  p  is  a  quantity  that  measures  the  strength  of 
the  interaction  of  the  collision  partner  with  the  diffuse  electron  cloud;  and 
y  is  a  quantity  related  to  the  energy  gap  between  the  nd  and  higher 
nA  states  of  Na.  The  functions  of  g  and  f  are  not  analytic  but  have  been 
calculated  and  make  simple  graphs.  This  equation  thus  leads  to  an  interpreta¬ 
tion  of  the  scattering  process  in  terms  of  reduced  parameters  that  have  a 
clear  physical  meaning. 

The  formula  has  been  applied  to  collisions  of  Na(nd)  with  He,  Ne,  Ar,  N2, 
CH^,  C^Hgj  to  Rb(nf)  +  He,  Ar,  Xe;  and  to  Xe(nf)  +  C02»  The  results  are 
generally  accurate  to  about  a  factor  of  two.  This  degree  of  accuracy  is  very 
useful  because  the  cross  sections  vary  by  an  order  of  magnitude  or  more  in  the 
range  of  n  considered.  A  paper  on  this  work  has  been  published  in  Physical 
Review  A  (see  Section  III  and  Appendix  A). 
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ABSTRACTS  OF  PUBLICATIONS 


A  coupled-channel  study  of  electronic-vibrational- 
translational  energy  exchange  in  nonadiabatic  collisions 

C.  H.  Becker  and  R.  P.  Saxon 

Molecular  Physics  Laboratory,  SRI  International.  Menlo  Park.  California  94025 
(Received  15  June  1981;  accepted  13  July  1981| 

Fully  quantum  mechanical  coupled -channel  calculations  of  electronically  nonadiabatic  collisions  in  (he 
presence  of  many  vibrational  states  have  been  performed  over  a  range  of  collision  energies  from  0.27  to  43.5 
eV.  The  calculations  were  based  on  simplified  model  potentials  for  K  +  0;,  considering  the  doublet  potential 
curves  arising  from  the  ground  electronic  neutral  K(*S)  +  and  ionic  K*(‘S)  +  Of (*//,)  states,  and 

assuming  isotropic  interaction.  Comparison  of  calculated  scattering  properties  and  recent  experimental  results 
at  hyperthermai  energies  for  K  +  02  indicates  that  the  simplified  model  interaction  nevertheless  gives  a  good 
physical  description.  This  work  investigates  the  collision  dynamics,  while  providing  a  reference  for  other  more 
approximate  theoretical  approches.  Processes  leading  to  neutral  and  ionic  products  are  discussed  for  this 
system  where  many  vibronic  channels  are  strongly  coupled. 

J.  Chem.  Phys.  75,  4899  (1981). 


Comparison  of  ion  pair  formation  in  the  systems  Ar*+I 
and  K+l2 

A.  P.  Hickman  and  Keith  T.  Gillen 

Molecular  Physics  Laboratory.  SRI  International,  Menlo  Park.  California  94025 
(Received  1  April  1980;  accepted  30  June  1980) 

A  simple  model  that  has  been  used  extensively  by  Los  and  co-workers  to  treat  ion  pair  formation  in  collisions 
of  alkali  atoms  with  diatomic  molecules  is  extended  to  include  continuum  coupling  via  a  competing  Penning 
ionization  channel.  This  extended  model  is  then  used  to  calculate  the  differential  cross  sections  for  ion  pair 
formation  for  the  system  Ar’+I,  over  the  energy  range  28-15*  eV  and  to  compare  with  a  previous  treatment 
of  K  +  I,.  In  the  absence  of  significant  competition  from  continuum  processes,  Ar'  is  expected  to  behave 
in  a  manner  similar  to  K,  since  the  active  electron  is  an  unpaired  4j  electron  in  both  cases.  We  perform 
model  calculations  for  Ar'+I,  to  investigate  the  effects  of  varying  the  potential  curves  and  charge 
exchange  matrix  elements  and  of  including  a  continuum  coupling  function  riR)  Comparison  with 
previous  calculations  for  K  +  I,  suggests  increased  repulsion  on  the  Ar’-I,  surfaces  relative  to  those  of 
|4-I,.  The  competing  mechanisms  of  excitation  transfer  and  Penning  ionization  may  have  a  small  effect 
upon  the  ion  pair  angular  distributions. 

J.  Chem.  Phys.  73,  3672  (1980). 


Model  for  fast,  nonadiabatic  collisions  between  alkali 
atoms  and  diatomic  molecules 

A.  P.  Hickman 

Molecular  Physics  Laboratory.  SRI  International.  Menlo  Park,  California  94025 
(Received  26  June  1980;  accepted  17  July  1980) 

Equations  for  collisions  involving  two  potential  surfaces  are  presented  in  the  impact  parameter 
approximation.  In  this  approximation,  a  rectilinear  classical  trajectory  is  assumed  for  the  translational 
motion,  leading  to  a  time-dependent  Schroedinger's  equation  for  the  remaining  degrees  of  freedom.  Model 
potentials  are  considered  for  collisions  of  alkali  atoms  with  diatomic  molecules  that  lead  to  a  particularly 
simple  form  of  the  final  equations.  Using  the  Magnus  approximation,  these  equations  are  solved  for 
parameters  chosen  to  model  the  process  Cs  +  0:— Cs * -f  Of,  and  total  cross  sections  for  ion-pair  formation 
are  obtained  as  a  function  of  energy.  The  results  exhibit  oscillations  that  correspond  qualitatively  to  those 
seen  in  recent  measurements.  In  addition,  the  model  predicts  that  the  oscillations  will  become  less 
pronounced  as  the  initial  vibrational  level  of  O,  is  increased. 

J.  Chem.  Phys.  T2,  4413  (1980). 
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Approximate  scaling  formula  for  colllsional  angular-momentum  mixing  of  Rydberg  atoms 

A.  P.  Hickman 

Molecular  Physics  Laboratory,  SRI  International  Menlo  Park.  California  94025 
(Received  27  June  1980) 

An  approximate  scaling  formula  has  been  determined  that  permits  the  rapid  estimation  of  crow  sections  for 
angular-momentum -changing  collisions  of  Rydberg  atoms  with  a  variety  of  targets,  using  information  about  low- 
energy-electron  scattering  from  the  target  The  formula  is  obtained  by  fitting  the  results  of  coupled-channel  and 
Bom-approximation  calculations  to  functions  of  reduced  parameters.  Application  to  Nalnd)  +  He,  Ne,  Ar.  N„  CH„ 

C>H|,  Rb(«/)  4-  He,  Ar,  Xe;  and  Xe)n/1  +  CO,,  suggests  that  the  accuracy  is  about  a  factor  of  2. 


Phys.  Rev,  A  23,  87  (1981). 


A- 2 


APPENDIX  B 


QUANTUM  MECHANICAL  CALCULATIONS  OF  SYMMETRIC  MOLECULAR  CHARGE  EXCHANGE 
AND  VIBRATIONAL  EXCITATIONS  IN  02+  +  02  COLLISIONS 


C.  H.  Becker 

Molecular  Physics  Laboratory 
SRI  International 
Menlo  Park,  California  94025 


ABSTRACT 


A  scattering  formulation  is  presented  for  symmetric  molecular  charge 

exchange  treating  the  internal  vibrational  degrees  of  freedom  and  translation 

quantum  mechanically.  A  unitary  transformation  is  found  to  reduce  the  2n 

coupled  radial  equations  to  two  sets  of  n  coupled  equations.  The  formulation 
+  2  3  - 

is  applied  to  0  (X  n  ,  v)  +  0,(X  £  ,  v' )  at  1,  8,  and  36  eV  center-of-mass 

L  g  t  g 

collision  energies.  Convergence  with  respect  to  vibrational  basis  set  size  is 
found  to  be  slower  at  lower  energy.  Integral  cross  sections  and  transition 
probabilities  as  a  function  of  nuclear  orbital  angular  momentum  are 
presented.  The  integral  cross  sections  summed  over  final  levels  for  vibra¬ 
tional  excitation  and  deexcitation  and  for  charge  transfer  are  on  the  order  of 
2 

10  A  .  Comparisons  of  earlier  classical  path  calculations  with  the  present 
results  do  not  show  quantitative  agreement  at  these  energies. 


MP  81-198R 
03/05/82 
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I.  INTRODUCTION 


Symmetric  atomic  and  molecular  charge  transfer  are  important  in  the 
behavior  of  ionized  gases.  While  a  great  deal  has  been  learned  about  symmet¬ 
ric  atomic  charge  transfer,^  comparatively  little  work  has  been  done  on 
symmetric  molecular  charge  transfer.  In  the  molecular  case,  the  exchange  of 
vibrational-translational  energy  as  well  as  electron  exchange  may  be  quite 
efficient. 

The  first  work  that  explicitly  included  the  vibrational  degrees  of 

2 

freedom  for  molecular  charge  exchange  was  by  Bates  and  Reid.  Ihis  was  a 
classical  path  (also  known  as  impact  parameter)  type  scattering  approximation 
using  a  time-dependent  framework.  The  translational  motion  is  treated 
classically  and  the  internal  degrees  of  freedom  are  treated  quantum 
mechanically.  This  treatment  is  expected  to  be  reliable  in  the  100-1000  eV 
collision  energy  range.  The  approximation  surely  degrades  at  some  point  below 
this,  but  it  is  unclear  where. 

Recently,  Moran  et  al.  performed  calculations  studying  O2  +  O2 
collisions  between  0.5  and  8  eV  kinetic  energy  with  the  formulation  of 
Ref.  2.  Comparison  was  made  between  a  straight-line  trajectory  and  an 
effective  curved  trajectory  for  the  classical  path.  It  is  reasonable  to 
expect  that,  for  lower  energies,  a  curved  trajectory  approach  would  give  a 
better  description. 

This  paper  presents  a  formulation  for  the  symmetric  molecular  charge 

exchange  process  that  treats  all  degrees  of  freedom  quantum  mechanically.  A 

unitary  transformation  is  then  found  to  reduce  the  dimensionality  of  the 

problem.  The  formulation  is  then  applied  to  <>2+  +  O2  collisions  and  the 

3—5 

results  are  compared  with  classical  path  calculations. 
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il.  THEORY 


A  quantum  mechanical  coupled-channel  formulation  is  used  to  describe  the 
symmetric  molecular  charge  exchange  scattering.  We  assume  that  the  symmetric 

scattering  is  rotationally  elastic,  as  assumed  in  the  earlier  treatment  of 

o 

Bates  and  Reid.  Therefore,  the  rotational  eigenfunctions  are  dropped  from 

the  basis  set,  though  they  could  be  included  in  a  straight-forward  manner. 

Rotational  excitation  could  in  principle  be  treated  exactly  or  treated  within 

some  approximate  framework  if  anisotropic  potentials  are  known. ^ 

Consider  two  molecules,  labeled  A  and  B,  with  vibrational  coordinates  ra 

and  rb,  and  the  distance  of  separation  of  the  centers-of-mass  of  the  two 

molecules  given  by  R.  In  theory,  we  can  distinguish  two  ground  state 

+  2 

electronic  configurations:  when  molecule  A  is  C>2  (X  IT^,  v)  and  molecule  B 
3 

is  02(X  ,  v'),  the  “direct,"  or  D,  electronic  state  results,  and  when 

molecule  A  is  (^(X^E  ,  v' )  and  B  is  02+(X^IIg,  v)  the  "exchange,"  or  X, 
electronic  state  results.  Higher  electronic  states  of  02+  and  02  are  ignored; 
this  should  not  have  a  significant  effect  at  the  investigated  collision 
energies. 

The  total  wavefunction  is  taken  as  the  diabatic  expansion 


D  D 

?(r  ,r  ,R,p)  -  <|,  (p)E  $  (r  )x  ,  (r  )F  (R) 

g  q  V  V  &  V  *  D  V 


+  (|>X(p)E  $  (r  )x  ,(r  )FX(R) 

V  V  D  V*  A  V 


(i) 


D  X 

where  <J»  ’  (p)  are  the  diabatic  electronic  eigenfunctions  of  the  electronic 
coordinates  p,  and  and  Xy» (r)  are  vibrational  eigenfunctions  for 
02+  and  02,  respectively. 
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Combinations  of  C>2+  and  0£  vibrational  levels  are  denoted  by 
D  X 

v  =  (v,  v'),  and  F  ’  (R)  are  the  radial  wavef unctions  for  each  channel 
containing  the  scattering  information  (at  large  R). 

The  total  Hamiltonian  is  written 


H  -  Hel  +  T(R)  +  T(ra)  +  T(rfa)  (2) 

where  Hfi^  is  the  Hamiltonian  for  the  electronic  degrees  of  freedom,  and  T  is 
the  kinetic  energy  operator.  The  projection  of  Hfi^  on  the  electronic 
eigenfunctions  yields 


Hel  /  "  u°(rHD  +  tVra>  +  (3a) 

Hel/  ’  /<*>/  +  tV0(ra}  +  *  (3b) 


D  X 

Here  U  *  represents  the  inte molecular  interaction,  and  Vq,  V+  are  the 
neutral  and  ionic  diatomic  potential  functions. 

Substituting  Eqs.  (1,  2,  and  3)  into  the  time-independent  nonrelativistic 

_  ,  D  *  X  * 

SchrMdinger  equation,  multiplying  on  the  left  by  (4>  )  ,  or  (<J»  )  ,  and 

integrating  over  p  yields 


[ V°(R)  +  T(R)  +  T(r  )  +  V  _<r  )  +  T(r)  +  V  (r.  )  -  E] 
a  +,0  a  b  0,+  b  J 


D,X 


*  £♦  (r  Jx  ,(rw  )F  (R)  + 
v  v  a,b  v1  b,a  v 


(A) 


[V^R)  +  C(T<R)  +  T(ra)  +  VQ  +(r#)  +  T(rb)  +  V+  Q(rb)  -  E)] 


x  £*  (r.  )X  , (r  )  l,D  (R)  -  0 
v  v  b,a  v  a,b  v 
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where 


C  -  <4-°|4'X>  “  <<PX|«I»D>,  VX(R)  =  <4,D|UX(p,  R)|*,S  -  «t,X|UD(p,  R)  1 4»D>p* 

V°(R)  =  <4)° |U°(p,  R)|4,D>p  -  <4»X j UX(p ,  R) 1 4»X>p *  and  Eq.  (4)  implies  the  two 
equations  corresponding  to  the  first  or  second  sub-  or  superscripts.  In 
writing  Eq.  (4)  we  have  neglected  the  nonadiabatic  coupling  of  the  nuclear 
kinetic  energy  operators  acting  on  the  electronic  wavef unctions,  that  is  we 
have  invoked  a  Born-Oppenheimer  approximation. 

Noting  the  vibrational  SchrBdinger  equations 

(T(ra,b>  +  Vra,b>)  ♦v^a.b5  ’  W'a.b*  and 

(T<ray  +  V0(ra,b)^V<ra,b)  "  W^a,^’  writin«  ev  "  ev  +  V* 
using  the  indices  v  and  p,  and  closing  on  Eq.  (4)  with  <$y(ra)  |  <xyt  (r.  ) | 

yields 

[T(R)  +  V°(R)  +  ep  -  E]F°(R) 

(5) 

+  5[V*(R)  +  C(T(R)  +  ev  -  BJ^FJCE)  -  0 

and 


[VX(R)  +  C(T(R)  +  «  -  E)]F®(R) 


+  E[v^ (R)  +  T(R)  +  e  -  e]p  FX(R)  -  0 

VL  V  J  pV  V 


(6) 


where  P  ■  P  ■  p(v  ,  v'  )p(v  ,  v'  )  in  which  p(v,  v' )  ■  <$  (r)|x  ,(r)> 
pv  vp  p  v  v  p  v  1  Av 

is  a  vibrational  overlap  integral. 

Rewriting  Eq.  (6)  as 


E[T(R)  +  e  -  E]P  FX(R) 
VL  V  J  pv  V 


-  -  §  VD(R)P^vFX(R)  -  [vX(R)  +  C(T(R)  +  cp  -  E)]F°(R) 


(7) 
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A. .. 


1 


1 


and  substituting  Eq.  (7)  into  Eq.  (5)  and  rearranging  yields 


{(1  -  C2)[T(R)  +  +  V°(R)  -  E]  -  C[VX(R)  -  CVD(R)]}f“(R) 


+  §[V*(R)  -  CV°(R)]P  F*(R)  -  0. 


For  symmetric  systems,  we  use  the  relationship' 


o  _  X  D 

(l  -  c  )v  -  v  -  cv 


where  V(R)  is  the  governing  intermolecular  interaction  (to  be  discussed 
shortly).  Using  Eq.  (9),  Eq.  (8)  is  rewritten  as 


[l(R)  +  e  +  V°(R)  -  CV(R)  -  E]F°(R)  -  -  EV(R)P  F*(R)  .  (10) 

r  r  ’  r''  ' 


Analogously,  one  finds 


[T(R)  +  e  +  V°(R)  -  CV(R)  -  E]FX(R)  -  -  §V(R)P  F°(R)  .  (l] 


For  n  combinations  of  v  and  v* ,  Eqs.  (10)  and  (11)  constitute  a  set  of  2n 
coupled  equations — n  direct  and  n  exchange  channels. 

For  numerical  solution  of  Eqs.  (10)  and  (11)  a  simple  partial  wave  expan¬ 
sion®  is  taken  for  FD,X(R), 


p  1  k  R  X  p 

P 


(12) 


where  is  Che  Legendre  polynomial,  and  6  is  Che  scaCCering  angle.  More 

precisely,  one  should  consider  Che  nuclear  statistics.^  For  che  nuclear 

spin  is  zero,  Bose-EinsCein  scacistics  govern,  and  Che  expansion  Eq.  (12) 

should  be  wriCCen  wich  Che  summation  2  £  . . .  .  For  large  numbers  of  parcial 

Jt(even) 

waves,  as  in  Che  presenC  application,  lictle  error  is  incurred  in  Che  integral 
cross  sections  using  Eq.  (12).  The  set  of  equations  of  the  form  of  Eqs.  (10) 
and  (11)  can  Chen  be  written,  now  in  matrix  form,  as 

2 

[I  ~  +  k2  -  I  ,  W(R)]f*(R)  .  q  .  (13) 

~  dR  ~  R 


This  is  a  matrix  equation  of  2n  dimension.  Matrix  elements  are  Indexed  by  p 

and  v  and  by  D  and  X.  Here  I  is  the  unit  matrix,  l  is  the  nuclear 

2 

orbital  angular  momentum,  k  is  a  diagonal  matrix  with  elements 
2  2 

k  -  — £  (E  -  e  )  6  and  W(R)  is  a  symmetric  matrix  with  elements 

pv  jf2  p  p,v  D,X 


DX 


V°  '  ^  1“  -  5D,X,P„v'(E)  +  6,,v6D,xivD(R)  -  •  <14> 

DX  * 


The  radial  wavefunctions  are  given  by  the  2n  dimensional  square  matrix 
f*(R);  each  column  is  a  linear  independent  solution  to  the  scattering 
problem  corresponding  to  a  different  initial  channel. 

I 

The  matrix  of  regular  solutions  f  (R)  can  be  written  as 


f  ■  u  - 


u+  k -1/2  S  k1/2 


(15) 


where  u*  are  diagonal  with  the  elements  being  the  Ricattl-Hankel  functions 

at  large  R,1®  and  approaching  exp[±  i  (k  R  -■—£)]  as  R  +  ®.  Also 

±1/2  ±1/2 

f(0)  -  0.  The  k  are  diagonal  with  elements  k  dropping  the 


B-7 


subscripts,  and  at  large  R  S  is  the  symmetric  and  unitary  scattering 
matrix,  ttie  transition  matrix  is  given  by  T  ■  I  -  S. 

IM  »V 

Using  Eq.  (15),  the  matrix  equation  (13)  can  be  written  schematically  in 
0  and  X  blocks  as 


ul0\  /u1o\/k'1/2lo  \/sddIsdx\A1/2I  0 

efW ; ! :-)  !?)  (t 


o  (16) 


Here  a,  u*,  and  k*^2  are  n-dimensional  diagonal  matrices,  a  has 
2 

elements  — s-  +  k2  -  [VD(R)  -  £V(R)],  b  has  symmetric 

ab2  W»  D2  *2  1  1 


SXD  S  SX*  and  SDD 


S  =  S  . 
XX  D 


R"  4i 

elements  P  V(R), 

li2 

The  symmetry  of  the  system  can  be  used  to  reduce  the  problem  from  2a 
dimensions  to  two  problems  each  of  n  dimensions.  As  the  computing  time  for 
numerical  solution  of  these  equations  generally  scales  as  the  total  number  of 
channels  cubed,  this  can  make  an  inqportant  difference.  Also  computer  storage 
problems  are  diminished.  The  unitary  matrix  which  diagonalizes  Eq.  (16)  is 
found  to  be  a  general  rotation  matrix,  written  in  block  form  as 


/ cos  (x/4)I  I  sin  (x/4)l\ 
(~sTn—  ( nj AjI  j  coJ  Tx7^Tl J 


\ 


y 


z 


(17) 


Multiplying  Eq.  (16)  on  the  left  by  ZC  (the  t  signifies  transpose)  and  on 
the  right  by  Z  leads  to  a  block  diagonal  form.  The  two  blocks  are  the  matlx 
equations 

(l±b){u--u+k"1/2(S<)±Sx)k1/2}  -0  .  (18) 
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By  integrating  Eq.  (18)  to  large  R,  for  both  signs,  the  Sp  and  Sx  then  can 
be  constructed. 

In  analogy  to  atom-atom  charge  (or  excitation)  transfer, ^  the  plus  solu¬ 
tion  of  Eq.  (18)  corresponds  to  gerade  (g)  state  scattering  while  the  minus 
solution  corresponds  to  ungerade  (u)  state  scattering.  The  S  and  S 

Is  A 

matrices  are  reconstructed  from  the  g  and  u  solutions. 

The  interaction  potential  V(R)  is  taken  as  one-half  the  gerade-ungerade 
eigenvalue  splitting.8  The  identification  of  V(R)  from  Eq.  (9)  as  one-half 
the  g-u  splitting  is  based  on  a  molecular  analogy  to  the  equivalence  between 
the  two  atom  LCAO  approximation  for  2V(R)  and  the  g— u  splitting,  neglecting 

Q 

the  electron  exchange  and  electron  motion  between  atoms.  For  these 
calculations  the  isotropic  V(R)  is  taken  from  Ref.  4. 

The  diagonal  potential  term  V0  -  £V  ia  found  from  the  Morse  and  Sato 

A 

function  evaluation  of  the  gerade  and  ungerade  eigenvalues 

eg  -  (V°  -  CV)  +  V  -  De{exp[-2p(R-Re)]-2exp[-P(R-Re)]}  (l9a) 

eu  “  (V0  "  CV)  -  V  -  |  De{exp[-2(J(R-Re)]  +2exp[-p(R-Re)]}  (19b) 

From  Eq.  (19)  we  find 

V°  ~  CV  -  D  {(3/4)exp[-2p(R-R  )]  -( l/2)exp[-p(R-R  )]}  .  (20) 

Due  to  the  Sato  approximation  for  eu»  there  is  an  unphyslcal  crossing  of 
eg  and  eu  at  R  (“4.306  bohr  for  the  parameters  used^).  For  values  of 
R<R  ,  V^-£V  is  taken  to  be  e  ,  and  V  is  taken  to  be  zero. 

s 


III.  RESULTS  AND  DISCUSSION 


Numerical  Integration  of  the  coupled  equations  (Eqs.  (18))  was  performed 
for  02+  +  O2  by  the  log-derivative  integration  algorithm  of  Johnson.**  For 
the  confutations  at  8  and  36  eV,  where  all  channels  of  the  basis  set  are  ener¬ 


getically  open,  the  integration  was  performed  out  to  a  distance  beyond  all 

classical  turning  points  where  a  connection  was  made  to  a  multichannel  JWKB 
1 2 

propagation  (see  also  Ref.  13).  By  propagating  the  multichannel  JWKB 

1  2 

solution  over  small  steps  using  the  repeated  Magnus  approximation,  the 
S-matrix  elements  from  the  JWKB  solution  converged  to  the  S-matrix  elements 


from  the  fully  integrated  equations,  at  least  within  the  working  accuracy  of 


the  integration  which  is  |S^|  values  accurate  to  within  0.001.  (The  log- 


derivative  integration  and  JWKB  propagation  step  sizes  and  end  points  were 

2 

chosen  to  maintain  an  accuracy  in  ls^vl  0.001  for  a  given  basis  set.) 

The  effects  of  different  basis  set  sizes  on  the  results  will  be  discussed 


shortly. 


The  products  of  overlap  integrals  were  calculated  by  direct  numerical 

integration  of  Morse  function  vibrational  wavef unctions.  The  Morse  parameters 

were  taken  from  Ref.  14  for  0„+(X^II  )  and  0„(x\  ).  With  regard  to 

2  g  2  g 

3—5 

comparisons  to  the  scattering  results  of  Moran  et  al. ,  the  overlap  integral 
products  used  in  this  work  do  not  differ  from  the  values  derived  from  RKR 


potentials  reported  in  Ref.  4  by  more  than  about  0.001.  This  difference  is 
expected  to  be  unimportant  in  comparisons  between  the  different  scattering 
computational  approaches. 
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Calculations  at  1  eV 

In  constructing  the  basis  set,  the  number  of  channels  was  incremented 
according  to  the  maximum  number  of  vibrational  quanta  in  the  0 ^  +  CU 

system.  For  example,  for  a  maximum  of  2  quanta,  the  number  of  channels  in  the 
basis  set  (ignoring  presently  0  versus  X  in  the  channel  count),  is  n  •  6 
corresponding  to  the  combinations  v  “  (v,  v' )  ■  (0,0),  (0,1),  (1,0),  (0,2), 
(1,1),  (2,0).  For  3  quanta  n  -  10,  for  11  quanta  n  *  78,  i.e.  ,  n  is  the 
binomial  coefficient  where  y  is  the  maximum  number  of  quanta. 

Scattering  calculations  at  1  eV  collision  energy  relative  to  0£+(v  ■  1)  + 
02(v*  -  0)  were  performed  at  only  one  value  of  orbital  angular  momentum 
(i.  «  92,  corresponding  to  an  impact  parameter  b  -  2.00  bohr)  because  it  was 
found  that  the  number  of  channels  required  to  achieve  well  converged  results 
precluded  calculations  over  many  values  of  JL.  A  study  of  the  vibrational 
basis  set  convergence  properties  at  1  eV  collision  energy  for  the  initial 
channel  (1,0)  is  presented  in  Table  l.  Values  of  the  square  modulus  of  the  D 
and  X  S-matrix  elements  are  presented.  Phase  convergence  of  the  S-matrix 
elements  follow  a  similar  trend.  The  slowest  channels  to  converge  are,  as 
expected,  the  elastic  D  and  X  channels,  while  the  energetically  closest 
channel  (0,1)  (AE  =  320  cm-^)  is  second-slowest  to  converge.  The  slowness  of 
the  elastic  channel  convergence  probably  reflects  the  sensitivity  of  the 
charge  transfer  process  to  modest  changes  in  scattering  phase  accumulation 
with  the  addition  of  more  coupled  channels.  As  can  be  seen  in  Table  1,  about 
78  channels  are  required  to  achieve  agreement  between  consecutive  basis  sets 
to  about  IX  for  the  elastic  D  and  X  probabilities. 

A  comparison  of  the  relative  probabilities  for  the  transitions  from  (1,0) 
for  b  -  2.00  bohr  and  1  eV  by  the  present  calculations  (Table  1)  with  the 
straight  line  and  curved  trajectory  classical  path  calculations  of  Ref.  3  show 
essentially  no  agreement  between  any  of  the  three  methods. 
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Calculations  at  8  eV 


Channel  convergence  tests  for  the  (1,0)  channel  for  kinetic  energy  of 

8  eV  relative  to  0£+(v  -  1)  +  02(v*  -  0)  were  performed  up  to  n  ■  45  channels 

(ignoring  D  versus  X  in  the  channel  count)  for  A  -  136  and  816.  It  was  found 

I  D  XI2 

that  n  »  21  channels  at  A  -  136  gave  convergence  of  S'  for  the 

I  H»v| 

slowest  converging  elastic  channels  to  within  about  4%,  and  within  1Z  for 
inelastic  channels,  which  was  deemed  sufficient  to  carry  out  the  calculation 
over  a  grid  of  A-waves  (inpact  parameters).  Better  convergence  was  found 
for  A  -  816.  The  spacing  used  was  AA  ■  34  corresponding  to 
Ab  -  0.26  bohr.  The  larger  number  of  channels  needed  at  the  lower  collision 
energy  perhaps  is  indicative  of  more  phase  accumulation  or  more  vibrationally 
adiabatic  behavior.  At  both  1  and  8  eV  the  O2"*",  O2  molecules  have  enough  time 

during  the  collision  for  more  than  one  vibrational  period.  Thus  classical 

13 

vibrational  timing  arguments  do  not  seem  capable  of  explaining  the 
convergence  difference. 

The  relative  transition  probabilities  as  a  function  of  A  or  b  are  seen 
plotted  in  Figs.  1  and  2  in  a  manner  that  allows  direct  comparison  to  the 
straight-line  and  effective  curved  trajectory  calculations  of  Ref.  3. 
Comparison  with  the  two  classical  path  approximations  shows  that  agreement 
with  the  fully  quantal  calculations  is  at  best  only  qualitative  for  this 
system  at  8eV. 

Absolute  transition  probabilities  for  the  scattering  at  8  eV  out  of  (1,0) 
are  shown  as  a  function  of  A  for  the  major  inelastic  direct  channels  in  Fig.  3 
and  inelastic  exchange  channels  in  Fig.  4.  The  elastic  exchange  probabilities 
are  plotted  in  Fig.  5.  Some  interesting  similarities  appear  between  several  D 
and  X  proabilities  in  Figs.  3  and  4,  but  a  sound  explanation  for  the  detailed 
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form  of  the  transitions  is  not  apparent  at  this  time.  Note  that  the 
(1,0)  +  (0,1)  transition  probability  extends  in  i.  considerably  beyond  the 
other  inelastic  channels  leading  to  larger  integral  cross  sections,  as 
expected.  Inelastic  transitions  as  a  function  of  i  involving  the 
(de)excitation  of  one  or  more  total  quanta  show  a  maximum  and  a  falloff  at 
l  *  500,  which  corresponds  to  a  classical  turning  point  around  R*,  the  e 

S 

and  ey  crossing  at  4.3  bohr.  That  the  critical  classical  turning  point  is 
around  R*  (this  is  also  found  at  36  eV),  and  not  greater  than  or  equal  to  the 
maximum  in  V  (5.0  bohr),  suggests  that  a  large  dV/dR  is  important  in 
producing  more  highly  inelastic  transitions,  because  for  R  just  greated  than 
R*  the  dV/dR  is  largest.  The  form  of  dV/dR  near  the  and  crossing  is 
strongly  model  dependent.  It  is  noted  that  the  cross  sections  for 
(1,0)  -*■  (1,0)  and  (0,1)  are  less  sensitive  to  this  region  of  e^,  eu  crossing. 

Values  of  some  integral  cross  sections  at  8  eV  from  02+(v  »  1) 

+  02(v  -  0)  are  given  in  Table  2  together  with  some  comparisons  to  classical 

3  Jl  2 

path  calculations.  Spline  interpolated  values  of  | T  |  are  used  in  the 

present  results,  so  the  values  are  approximate  but  should  be  reliable  to 

within  a  few  per  cent  as  the  figures  indicate.  The  figures  also  indicate  that 

the  chosen  spacing  in  A£  gives  an  accurate  description  of  the  collision 

because  of  the  smoothness  of  the  probabilities  as  a  function  of  1.  The 

results  of  Table  2  and  comparison  of  Figs.  1  and  2  with  those  of  Ref.  3  do  not 

indicate  that  the  curved  trajectory  is  a  significant  improvement  over  a 

straight-line  trajectory.  Agreement  is  not  particularly  close  between  any  of 

the  three  approaches  at  8  eV. 
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Calculations  at  36  eV 


Channel  convergence  tests  for  scattering  at  36  eV  center-of-mass  colli¬ 
sion  energy  relative  to  the  (0,0)  channel  showed  that  the  basis  set  conver¬ 
gence  was  more  rapid  than  at  8  eV.  A  total  of  n  ■  21  channels  were  used  In 
the  coupled -channel  calculations  with  a  grid  of  i  waves  using  M  •  68 
corresponding  to  Ab  *  0.23  bohr.  Absolute  transition  probabilities  as  a 
function  of  l  for  scattering  from  (1,0)  are  shown  for  the  elastic  charge 
exchange  process  in  Fig.  3,  for  the  inelastic  direct  channels  in  Fig.  6,  and 
for  the  inelastic  charge  exchange  channels  in  Fig.  7.  Comparison  of  the  8  and 
36  eV  elastic  charge  exchange  probabilities  (Fig.  3)  shows  significant 
similarities;  the  main  differences  are  that  transitions  occur  out  to  slightly 
larger  impact  parameters  at  lower  energy  and  that  the  8  eV  results  show 
additional  strong  oscillations  (phase  interference),  both  suggesting  the 
importance  of  more  phase  accumulation  at  lower  energy. 

Total,  i.e.,  summed  over  final  states,  inelastic  and  total  charge 
exchange  integral  cross  sections  Q  at  36  eV  from  the  present  results  and  those 
from  straight-line  trajectory  calculations  of  Refs.  4  and  5  are  compared  in 
Table  3.  As  can  be  seen,  agreement  is  qualitative  only.  Similar  agreement  is 
found  for  state-to-state  integral  cross  sections.  It  is  expected  that  as  the 
energy  increases,  the  classical  path  approximation  will  be  increasingly 
accurate. 

Charge -exchange  and  vibrationally  inelastic  cross  sections  have  been 
found  to  be  rather  large,  as  seen  in  Tables  2  and  3  indicating  rapid  energy 
migration  in  ionized  oxygen,  and  large  momentum  transfer  cross  sections. 

In  summary,  a  quantum  mechanical  coupled-channel  formulation  has  been 
presented  for  symmetric  molecular  charge  exchange.  A  unitary  transformation 
was  found  to  reduce  the  set  of  2n  coupled  radial  equations  to  two  sets  of  n 
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coupled  equations.  Calculations  have  been  applied  to  the  02+  +  02  system  at 
1,  8,  and  36  eV  center-of-mass  collision  energies.  Vibrational  basis  set  size 
convergence  is  very  slow  at  1  eV,  but  becomes  more  rapid  at  8  eV,  and  more 
rapid  still  at  36  eV.  Conparisons  of  classical  path  scattering  approximation 
calculations  with  the  present  calculations,  which  treat  translation  quantum 
mechanically,  show  substantial  quantitative  disagreement  for  the  02+  +  02 
system  at  these  energies.  Better  agreement  between  coupled-channel  and 
classical  path  calculations  would  be  expected  at  higher  collision  energies. 
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Table  1 


Values  of  Che  square  modulus  of  direct  and  exchange  S-matrix 
elements  for  selected  final  channels,  n,  from  the  v  -  (1,0) 
channel  at  1  eV  collision  energy  for  Z  -  92  as  a  function  of 
the  number  of  vibrational  channels  in  the  basis  set,  n. 


n 

|s°  I2 

|sx  |2 

|sD  |2 

1  (iV1 

|sx  |2 

'  (XV1 

n 

-  15 

n  -  21 

0,0 

0.00012 

0.00007 

0.00010 

0.00006 

0,1 

0.0081 

0.0930 

0.0127 

0.0861 

1,0 

0.1801 

0.7175 

0.1312 

0.7685 

0,2 

0.00080 

0.00009 

0.00076 

0.00007 

n 

-  28 

n  ■  36 

0,0 

0.00009 

0.00005 

0.00009 

0.00004 

0,1 

0.0174 

0.0690 

0.0207 

0.0558 

1,0 

0.1679 

0.7442 

0.2307 

0.6914 

0,2 

0.00075 

0.00006 

0.00072 

0.00005 

n 

-  45 

n  -  55 

0,0 

0.00008 

0.00003 

0.00008 

0.00003 

0,1 

0.0223 

0.0478 

0.0228 

0.0435 

1,0 

0.2886 

0.6400 

0.3288 

0.6035 

0,2 

0.00071 

0.00005 

0.00070 

0.00005 

n 

■  66 

n  -  78 

0,0 

0.00008 

0.00003 

0.00008 

0.00003 

0,1 

0.0229 

0.0414 

0.0229 

0.0404 

1,0 

0.3513 

0.5830 

0.3629 

0.5724 

0,2 

0.00069 

0.00005 

0.00069 

0.00005 
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Table  2  Integral  cross  sections  Q(A  )  at  8  eV  for 

02+(v  -  1)  +  02(v'  -  0)  -*•  02+(v)  +  02(vr)  Direct 

02(v')  +  02+(v)  Charge  Exchange 


Final  Channel  (v,  v' ) 

£. 

st 

0,0 

0.11,a 

1  0.66 

,b  0. 76c 

0.14, a 

1  0.37 

,b  o.: 

0,1 

2.23, 

1.75, 

1.79 

1.80, 

1.72, 

2.41 

1,0 

73.95 

26.40 

0,2 

0.90, 

0.55, 

0.64 

1.16, 

0.98, 

0.20 

1,1 

1.05, 

0.83, 

0.87 

1.23, 

0.93, 

1.52 

2,0 

0.03, 

0.48, 

0.70 

0.05, 

0.78, 

0.85 

0,3 

0.07, 

0.07 

0.04, 

0.05 

1,2 

0.02 

0.10 

2,1 

0.02 

0.02 

3,0 

0.01 

0.01 

a  Tlie  first  column  is  the  present  results. 

b  The  second  column  is  the  straight-line  trajectory  results  of  Ref.  3. 
c  The  third  column  is  the  curved  trajectory  results  of  Ref.  3. 


2 

Table  3  Integral  cross  section  Q(A  )  at  36  eV  for 
02+(v)  +  02(v')  -►  02+(v)  +  02(v')  Direct 

+■  02(v')  +  02+(v)  Charge  Exchange 


Initial  Channel  (v,  v' ) 

Total  Inelastic  Q 

Total  Charge  Exchange  Q 

0,0 

5.5 ,a  2.7b 

16. 4, a  29. b 

1,0 

13.6,  7.1 

29.5,  c 

2,0 

19.4,  11.1 

23.2,  31. 

aFlrst  column  Is  from  the  present  fully  quantum  mechanical  coupled-channel 
results. 

^Second  column  Is  from  straight-line  classical  path  calculations  of  Refs.  4 
and  5. 

cNot  available. 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 
Figure  3 


Figure  4 
Figure  5 


Figure  6 
Figure  7 


Relative  transition  probabilities  for  formation  of  specific  product 
channels  for  charge  transfer  (exchange)  scattering  at  8  eV  center-of- 
mass  energy  from  02+(v  *»  1)  +  02(v'  *  0).  Product  channels 
associated  with  each  area  are  denoted  by  (v,  v').  Calculations  by 
the  coupled-channel  method  were  performed  at  selected  values  of 
orbital  angular  momentum  1  given  by  the  symbol  +.  The  connecting 
lines  are  drawn  by  spline  Interpolation. 

Same  as  Figure  1,  except  for  direct  scattering. 

Absolute  values  of  direct  transition  probabilities  as  a  function  of 
£  out  of  (1,0)  at  8  eV  to  several  indicated  inelastic  channels. 
Calculated  values  are  marked  by  +.  The  connecting  lines  are  drawn  by 
spline  interpolation. 

Same  as  Figure  3  except  for  charge  exchange  scattering. 

Absolute  values  of  the  02+(v  »  1)  +  02(v'  -  0)  elastic  charge 
exchange  transition  probability  as  a  function  of  X  (and  b)  for  8 
and  36  eV.  Calculated  values  are  marked  by  +.  The  connecting  lines 
are  drawn  by  spline  interpolation. 

Same  as  Figure  3  except  for  36  eV. 

Same  as  Figure  3  except  for  36  eV  and  charge  exchange  scattering. 
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